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A coding system under consideration for the Galileo S-band low-gain antenna 
mission is a concatenated system using a variable redundancy Reed-Solomon outer 
code and a ( 14,1/4 ) convolutional inner code. The 8-bit Reed-Solomon symbols are 
interleaved to depth 8, and the eight 255-symbol codewords in each interleaved block 
have redundancies 64, 20, 20, 20, 64, 20, 20, and 20, respectively (or equivalently, the 
codewords have 191, 235, 235, 235, 191, 235, 235, and 235 8-bit information symbols, 
respectively). This concatenated code is to be decoded by an enhanced decoder 
that utilizes (1) a maximum likelihood (Viterbi) convolutional decoder; (2) a Reed- 
Solomon decoder capable of processing erasures; (3) an algorithm for declaring 
erasures in undecoded codewords based on known erroneous symbols in neighboring 
decodable words; (4) a second Viterbi decoding operation (redecoding) constrained 
to follow only paths consistent with the known symbols from previously decodable 
Reed-Solomon codewords; and (5) a second Reed-Solomon decoding operation using 
the output from the Viterbi redecoder and additional erasure declarations to the 
extent possible . 

It is estimated that this code and decoder can achieve a decoded bit error rate 
of 1 x 10~ 7 at a concatenated code signal- to-noise ratio of 0.76 dB. By comparison, 
a threshold of 1.17 dB is required for a baseline coding system consisting of the 
same ( 14,1/4 ) convolutional code, a (255,223) Reed-Solomon code with constant 
redundancy 32 also interleaved to depth 8, a one-pass Viterbi decoder, and a Reed- 
Solomon decoder incapable of declaring or utilizing erasures. The relative gain of 
the enhanced system is thus 0.41 dB. 

It is predicted from analysis based on an assumption of infinite interleaving that 
the coding gain could be further improved by approximately 0.2 dB if four stages 
of Viterbi decoding and four levels of Reed-Solomon redundancy are permitted. 
Confirmation of this effect and specification of the optimum four-level redundancy 
profile for depth-8 interleaving is currently being done. 
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I. Introduction 

This article looks at the performance of a Reed-Solomon 
plus convolutional concatenated coding system with en- 
hanced decoding as planned for the Galileo S-band low- 
gain antenna (LG A) mission. The baseline system with- 
out enhanced decoding uses a (255,223) Reed-Solomon 
outer code concatenated with a (14,1/4) convolutional in- 
ner code, and interleaves the Reed-Solomon symbols to 
depth 8. The convolution ally encoded symbols are de- 
coded by maximum likelihood (Viterbi) decoding, and 
each Reed-Solomon codeword is decoded algebraically to 
correct a maximum of sixteen 8-bit symbols per 255- 
symbol word. There are two types of decoding enhance- 
ments planned for the Galileo LGA mission: Reed-Solomon 
redecoding using erasure declarations and Viterbi redecod- 
ing using Reed-Solomon corrected symbols. 

Reed-Solomon redecoding is possible when at least one 
but fewer than eight of the codewords within a block of 
eight interleaved words is decodable (correctable). The 
Reed-Solomon decoder can then extrapolate the locations 
of corrected errors in the decodable word(s) to neighbor- 
ing locations in adjacent undecodable word(s) and declare 
the corresponding symbols to be erased. The basis for this 
error extrapolation or error forecasting is that errors from 
the Viterbi decoder tend to occur in bursts that are of- 
ten longer than the 8-bit Reed-Solomon symbols. If the 
erased symbols are highly likely to be erroneous, then the 
undecoded words might be decoded by a second try at 
Reed-Solomon decoding that utilizes the erasure informa- 
tion. 

Viterbi redecoding consists of an extra pass through 
a maximum likelihood decoder that is now constrained 
to follow only paths consistent with the known symbols 
from previously decodable Reed-Solomon codewords. The 
Viterbi redecoder is much less likely to choose a long erro- 
neous path because any path under consideration is pinned 
to coincide with the correct path at the location(s) of the 
known symbols. 

Both redecoding processes may be repeated an arbi- 
trary number of times (with diminishing returns). Rep- 
etitions of Reed-Solomon redecoding allow the testing of 
several different combinations of reasonable error extrap- 
olations from the known error locations, and can be ac- 
complished relatively cheaply, as Reed-Solomon decoding 
is much faster than Viterbi decoding. Each repetition of 
Viterbi redecoding requires a whole new decoding stage: 
it must begin with the output from the previous Reed- 
Solomon decoding stage and must feed its output to an- 
other following Reed-Solomon stage. 


With both types of redecoding, it usually pays to put 
different amounts of redundancy in neighboring Reed- 
Solomon codewords. Words with high redundancy can be 
counted on to decode during an initial decoding try, and 
the information from these decoded words can be used to 
assist the decoding of codewords with lower redundancy 
later. 

In this article, the stages of the enhanced decoding pro- 
cess are denoted as follows: 

VIT-1 the first pass of the raw data through the Viterbi 
decoder. 

RS-1 the first pass of the Viterbi decoded bits through 
the Reed-Solomon decoder, including the possi- 
bility of several trials per codeword using era- 
sure declarations forecast from previously de- 
coded codewords. 

VIT-2 the second pass of the raw data through the 
Viterbi decoder, aided this time by known bits 
from codewords successfully decoded in RS-1. 

RS-2 the pass of the Viterbi redecoded bits from 
VIT-2 through the Reed-Solomon decoder, again 
using erasure declarations to assist in decoding 
previously undecodable codewords. 


VIT-n the nth pass of the raw data through the Viterbi 
decoder, aided this time by known bits from 
codewords successfully decoded in RS-(n — 1). 

RS-n the pass of the Viterbi redecoded bits from 
VIT-n through the Reed-Solomon decoder, 
again using erasure declarations to assist in de- 
coding previously undecodable codewords. 


II. Three Analysis Approaches 

The analysis in this article is first performed using an 
assumption of infinite interleaving and no Reed-Solomon 
erasure declarations. This is done for one, two, and four 
Viterbi decoding stages. A calculation is made of the max- 
imum coding gain obtainable by adding the extra stages. 
With the infinite interleaving assumption, concatenated 
system performance can be accurately obtained to error 
rates of 1 x 10” 7 or lower, based on megabits of simulated 
Viterbi decoder data. 
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The analysis then continues for the actual Galileo con- 
ditions with depth-8 interleaving but still no Reed- 
Solomon erasure declarations. This analysis cannot di- 
rectly verify 1 x 10“ 1 * * * * * 7 error rates even with gigabits of 
Viterbi decoded data, but extrapolated performance 
curves for depth-8 interleaving can be relied on with high 
confidence because of their similarity to certain corre- 
sponding curves derived under the infinite interleaving as- 
sumption. 

Finally, the depth-8 interleaving analysis is extended 
to allow Reed-Solomon erasure declarations using a simple 
error forecasting strategy. 

The analysis in this article characterizes the concate- 
nated decoding system based on the usual assumption that 
the received symbols are corrupted by stationary additive 
white Gaussian noise. No attempt has been made to re- 
fine the decoder models to account for effects due to such 
factors as imperfect arraying, phase drifts, or time- varying 
signal-to-noise ratio. 


III. Analysis Based on Infinite Interleaving 
and No Erasure Declarations 

This section analyzes the gains possible from Viterbi 
redecoding alone, without Reed-Solomon erasure declara- 
tions, but using an assumption of infinite interleaving. It 
begins with a baseline analysis of one-stage decoding (no 
redecoding) and proceeds to evaluate the gains possible 
from using two stages and four stages. 

A. One-Stage Decoding 

One-stage decoding for this case consists of one pass 
through the Viterbi decoder followed by one pass without 
erasure declarations through the Reed-Solomon decoder. 
The performance of the Viterbi decoder is evaluated by 
simulation, and that of the Reed-Solomon decoder is cal- 
culated from a formula. 

1. Stage VIT-1. The Viterbi decoder for the Galileo 

(14,1/4) convolutional code is characterized by the perfor- 

mance curves given in Fig. 1. The curves show Viterbi 

decoder bit error rate BER and 8-bit Reed-Solomon sym- 
bol error rate SER as functions of the signal-to-noise ra- 

tio E b /N 0 of the convolutional code. The points on these 

curves are based on about 2 Gbits of data each in the range 
from 0.0 to 0.5 dB, and on more than 400 Mbits each in 
the range from 0.6 dB to 1.0 dB (but 2 Gbits at 0.7 dB). 
All of the data were decoded by the Big Viterbi Decoder 

(BVD). 


Figure 1 shows that the (14,1/4) Galileo LG A code 
achieves a BER of 5 x 10 -3 at an E b /N 0 of 0.65 dB. By 
comparison, Tables 1 and 2 of [1] show that the (15,1/4) 
Galileo high-gain antenna (HGA) code achieves a BER of 
5 x 10“ 3 at an E b /N 0 of 0.52 dB. The (14,1/4) LGA code 
is 0.13 dB inferior to the (15,1/4) HGA code at this error 
rate. 

2. Stage RS-1. Figure 2 shows the performance of 
the Reed-Solomon decoder under an assumption of infinite 
interleaving and no erasure declarations, but with varying 
amounts of error correction E (or codeword redundancy 
2 E). The E values for the various curves are labeled along 
the right and bottom edges of the graph. These curves 
are based on the same BVD simulation runs described for 
Stage VIT-1. On this figure, and on all succeeding figures, 
the signal-to-noise ratio E b /N 0 on the x-axis is the same 
as that in Fig. 1, i.e., the signal-to-noise ratio for the inner 
convolutional code alone. Thus, these curves cannot be 
used for a direct reading of the optimum value of E for 
single-stage Reed-Solomon decoding because they do not 
reflect the increasing dilution of signal energy due to in- 
creasing redundancy of the Reed-Solomon code. However, 
this latter adjustment is a simple numerical computation, 
and presenting a parametric family of performance curves 
as those in Fig. 2 leads to effective trial-and-error opti- 
mization methods for the multistage decoding processes 
analyzed below. 

The curves in Fig. 2, and in all succeeding figures, show 
the Reed-Solomon decoded 8-bit SER rather than the de- 
coded BER. This is a convenience because the analysis of 
future decoding stages depends primarily on SER rather 
than BER . The BER is slightly less than half the corre- 
sponding SER, as can be seen by comparison of the BER 
and SER curves in Fig. 1. 

3. Performance of One-Stage Decoding. For the 
baseline coding system using the (255,223) Reed-Solomon 
code, E is 16 and a one-stage decoding system with infinite 
interleaving can achieve an SER of 2x 10" 7 (i.e, a BER of 
~ 1 x 10“ 7 ) at convolutional code E b /N 0 = 0.52 dB. The 
corresponding concatenated code signal-to-noise ratio is 
E b /N 0 — 1.10 dB, accounting for the 0.58 dB redundancy 
of the (255,223) Reed-Solomon code. 

Alternatively, if a (255,231) Reed-Solomon code with 
E = 12 had been used, the required convolutional code 
operating point would move to E b /N 0 = 0.67 dB, but 
the corresponding concatenated code signal-to-noise ratio 
would stay at E b /N 0 = 110 dB, because the redundancy 
of the (255,231) Reed-Solomon code is only 0.43 dB. In 
fact, the lowest required concatenated code signal-to-noise 
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ratio of approximately 1.10 dB is achieved over a range of 
values of E from 12 to 16, or equivalently over a range of 
convolutional code operating points Eb/N 0 from 0.52 dB 
to 0.67 dB. 

B. Two-Stage Decoding 

Two-stage decoding consists of two passes each through 
the Viterbi decoder and the Reed-Solomon decoder. The 
only difference in the decoding algorithms is that the 
Viterbi decoder on the second pass is constrained to fol- 
low only paths consistent with known symbols from Reed- 
Solomon codewords that were decoded on the first pass. 
The Reed-Solomon decoding algorithm is unchanged for 
both passes. 

1. Stage VIT-1. The output of Stage VIT-1 is char- 
acterized by the same performance curves given in Fig. 1. 

2. Stage RS-1. The output of Stage RS-1 is char- 
acterized by the same performance curves given in Fig. 2. 
However, now it is important to utilize the full set of in- 
formation contained in all the curves in Fig. 2 rather than 
just to focus on the values of E from 12 to 16 that give 
optimum one-stage performance. 

3. Stage VIT-2. The performance of Stage VIT-2 is 
dependent on how often known 8-bit symbols from the pre- 
vious Reed-Solomon decoding are inserted into the Viterbi 
decoder data stream. If the Reed-Solomon code has con- 
stant redundancy, the known symbols will occur with a 
certain average frequency (determined by the correspond- 
ing SER curve in Fig. 2), but not according to any repeat- 
ing pattern. If the Viterbi redecoder had to contend with 
only randomly occurring unknown symbols interspersed 
with probability SER < 10~ 1 into a stream of mostly 
known symbols, Stage VIT-2 would produce an almost 
perfect output at any Eb/N 0 shown in Fig. 2. However, 
even with infinite interleaving, Stage RS-1 decoding fail- 
ures will be correlated from one codeword to the next due 
to the burstiness of the errors from Stage VIT-1. There- 
fore, the Viterbi redecoder is likely to encounter strings of 
unknown symbols of unpredictable length. Long strings of 
unknown symbols will locally negate the value of the oth- 
erwise highly likely known symbols and cause correspond- 
ingly long strings of redecoded errors. No effective ana- 
lytical techniques have been developed to quantitatively 
assess the performance of Stage VIT-2 after decoding a 
constant redundancy code in Stage RS-1. 

If the Reed-Solomon codewords have variable redun- 
dancy and there is an appreciable difference in redun- 
dancy values, then it is possible (and reasonable) to set 


the Eb/N 0 operating point so that the most highly redun- 
dant codewords are almost certain to decode, while the 
less redundant codewords are undecodable with probabil- 
ity much higher than the target error probability for the 
multistage decoding process. Then, to a good conservative 
approximation, it can be assumed that the known symbols 
from Stage RS-1 will occur in at least the same repeat- 
ing pattern as the pattern of the most highly redundant 
codewords. Under these conditions, the performance of 
Stage VIT-2 can be effectively parameterized by feeding 
the Viterbi redecoder various repeating patterns of known 
symbols. 

Viterbi redecoded BER and SER tables were pre- 
sented in [2] for three different known symbol patterns 
(one in two, one in four, and one in eight) and three differ- 
ent convolutional codes (the (7,1/2) NASA standard code, 
an (11,1/4) code, and the (15,1/4) Galileo HGA code). A 
similar analysis for the (14,1/4) Galileo LG A code was 
performed for this article, and the SER results are shown 
in Fig. 3. These curves are based on slow-running soft- 
ware simulations of only about 1 Mbit each, because the 
high-rate BVD used to generate the 2-Gbit samples for 
Stage VIT-1 could not without modification incorporate 
the known symbol information into its decoding process. 

Figure 3 shows multiple curves for the one-in- four and 
one-in-eight known symbol cases, because, as noted in [2], 
the redecoded SER depends on the symbol’s phase relative 
to the position of the nearest known symbol in the pattern. 
It can be argued theoretically and has been confirmed by 
simulation that phases ±n have approximately the same 
SER ) so the curves in Fig. 3 represent averages of the 
simulated SER values for ±n (n = 1, 2 for the one-in-four 
case and n = 1, 2, 3, 4 for the one-in-eight case). 

Figure 3 has two different labels for its y-axis. The label 
on the left is the actual SER observed from the simulation 
runs. The label on the right is the equivalent convolutional 
code EbfN 0 that produces the same SER values without 
any known symbols, i.e., the convolutional code signal-to- 
noise ratio that would have been required to produce the 
same SER with one-stage Viterbi decoding. 

4. Stage RS-2. The output of Stage RS-2 is char- 
acterized by the same family of performance curves given 
in Fig. 2. Under the infinite interleaving assumption, all 
symbols in a given codeword are independent, so the de- 
coding error rate depends only on the SER from Viterbi 
decoding during Stage VIT-2. Thus, the curves in Fig. 2 
are directly applicable to the determination of the perfor- 
mance of Stage RS-2 if the equivalent one-stage Eb/N 0 
values are used to characterize the output of Stage VIT-2. 
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For example, suppose that every fourth codeword is 
strong (i.e., highly redundant) and can be presumed to 
decode with high probability during Stage RS-1. Then for 
Stage VIT-2, there are known symbols spaced four sym- 
bols apart, and the curves labeled 1/4 and 2/4 in Fig. 3 
are applicable. If the convolutional code operating point is 
E b /N 0 = 0.20 dB, then the equivalent Stage VIT-1 operat- 
ing point is 0.85 dB for symbols at phases dbl and 0.77 dB 
for symbols at phase i2. If a Reed-Solomon code with 
E — 9 is used for the remaining weak (i.e., less redundant) 
codewords, Fig. 2 shows that the SER for the weak code- 
words is approximately 1 x 10~ 7 for codewords at phases 
=bl and 1 X 10“ 6 for codewords at phase ±2. 

5. Performance of Two-Stage Decoding. The 
overall performance of two-stage decoding depends on the 
likelihood of decoding both the strong codewords during 
Stage US-1 and the weak codewords during Stage RS-2... 
The SER of the strong codewords can be conservatively 
presumed to equal the SER from Stage RS-1 for a code 
with the redundancy of the strong codewords. There are 
two contributions to the SER of the weak codewords. 
Weak codewords will fail to decode when the number of 
symbol errors from the re decoding stage exceeds the weak 
codewords’ correction capability. Furthermore, weak code- 
words may not successfully decode whenever neighboring 
strong codewords do not decode during Stage RS-1. In 
this article, it is assumed that this second contribution to 
the SER of weak codewords equals the SER of the strong 
codewords. This leads to an approximate formula for the 
overall SER of the form SER = SER s { 1) + f w SER w ( 2), 
where SER s (\) is the strong codewords’ symbol error rate 
after Stage RS-1; SER W (2) is the weak codewords’ symbol 
error rate after Stage RS-2, assuming known symbols dur- 
ing Stage VIT-2 from neighboring decodable strong code- 
words; and f w is the fraction of weak codewords. If the 
weak codewords do not all have equivalent performance, 
the second term of this formula can be generalized to in- 
clude a sum of contributions from each individual type of 
weak codeword, weighted by the fraction of each type. 

If the target SER is 2 x 10~ 7 , then both SER s (l) and 
SER W { 2) should be reduced to roughly 1 x 10~ 7 . Build- 
ing upon the example introduced in the previous section 
with one strong codeword every four, Fig. 2 shows that 
a Stage RS-1 correction capability of E = 30 is required 
to bring SER s ( 1) to 1 x 10“ 7 at a convolutional code op- 
erating point of E b /N 0 = 0.20 dB. Similarly, the discus- 
sion of the previous example noted that SER W ( 2) could 
be brought to 1 x 10“ 7 for the weak codewords at phases 
±1 by using a Reed-Solomon code with a correction ca- 
pability of E — 9. For the weak codewords at phase ±2, 
E = 9 is insufficient because the resulting SER W (2) of 


1 x 10“ 6 would raise the overall SER above 2 x 10 -7 all 
by itself. However, E = 10 yields SER W ( 2) = 2 x 10“ 7 
for the weak codewords at phase ±2, resulting in an over- 
all SER = 1 x 10- 7 + (1/2)1 x 1(T 7 + (1/4)2 x 1(T 7 = 

2 x 10 -7 . The redundancy profile for this codeword set 
is 2 E = (60,18,20,18), which corresponds to an average 
code rate of 229/255, or 0.52 dB of overhead for the Reed- 
Solomon code. The corresponding concatenated code op- 
erating point is E b /N 0 = 0.20 dB + 0.52 dB = 0.72 dB. 
This represents an improvement of 0.38 dB relative to the 
required operating point for the one-stage decoder. 

If the target SER is still 2 x 10" 7 but the convolu- 
tional code operating point is moved to 0.10 dB, the re- 
quired strong codeword correction capability to achieve 
SER = 1 x 10' 7 in Stage RS-1 is now E - 36. Af- 
ter Viterbi redecoding with one in four symbols known, 
the SE ^-equivalent Stage VIT-1 operating points from 
Fig. 2 are 0.79 dB for phases ±1 and 0.72 dB for phase 
±2. Weak-codeword correction capabilities of E = 10 and 
E = 11 are sufficient to bring the corresponding SER W ( 2) 
values to 1 x 10" 7 and 2 x 10" 7 , respectively. The overall 
SER = lx 10 -7 +(l/2)l x 10~ 7 +(l/4)2x 10 -7 = 2xl0~ 7 , 
as before. The redundancy profile of this codeword set 
is 2 E = (72,20,22,20), which corresponds to an average 
code rate of 443/510, or 0.61 dB of overhead for the Reed- 
Solomon code. The corresponding concatenated code op- 
erating point is E b /N 0 — 0.10 dB 4- 0.61 dB = 0.71 dB, 
which is almost exactly the same as before. As for the case 
of one-stage decoding, there is a range of convolutional 
code operating points over which the same near-optimal 
concatenated code performance can be achieved by appro- 
priately adjusting the Reed-Solomon redundancies. 

The case of one strong codeword every two is not di- 
rectly analyzable from Fig. 2 because the equivalent con- 
volutional code Eh /N 0 is outside the range of values mea- 
sured for Fig. 2. However, calculations show that for SER 
in the range of 0.0006 to 0.0009, a Reed-Solomon code with 
E = 4 will bring the SER for Stage RS-2 under 1 x 10~ 7 , 
but with E = 3 it is above 3 x 10~ 7 . Referring to Fig. 3 for 
the Stage VIT-2 SER values after redecoding with one in 
two symbols known, E = 4 will be required for the weak 
codewords if the convolutional code operating point is be- 
yond about 0.1 dB. Because of the large spacing between 
the E — 4 and the E = 3 curves, it can be argued that 
E = 3 remains inadequate for a few tenths of a dB beyond 
0.1 dB (and beyond the region in Fig. 3 for which simu- 
lation results have been obtained thus far). The choice of 
E for the strong codewords depends on the convolutional 
code operating point. For example, E = 30 is required 
at 0.20 dB, while E = 25 is sufficient at 0.30 dB. The 
corresponding redundancy profiles (60,8) and (50,8) con- 
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tribute Reed-Solomon coding overheads of 0.62 dB and 

0.52 dB, respectively, and both of the resulting concate- 
nated code operating points are E b /N 0 = 0.82 dB. This 
is about 0.1 dB worse than the best operating point ob- 
tainable with the scheme that uses only one strong code- 
word every four. It reflects the inefficiency of investing 
high amounts of Reed-Solomon redundancy in codewords 
spaced so closely together that the Stage VIT-2 SER in 
Fig. 3 is driven lower than necessary. 

The case of one strong codeword every eight blurs the 
distinction between strong and weak codewords because of 
the relatively wide variation in Stage VIT-2 SER among 
the various possible symbol phases. For example, at 
0.20 dB a code with E = 30 will bring the strong codeword 
SER to 1 x 10“ 7 . From Fig. 3, the SER-e quivalent oper- 
ating points for the weak codewords after Viterbi redecod- 
ing are 0.67 dB, 0.51 dB, 0.43 dB, and 0.41 dB for phases 
±1, db2, ±3, and ±4, respectively, requiring E = 12, 17, 19, 
and 20, respectively, to bring the individual SER W { 2) val- 
ues to the range of 1 x 10" 7 to 2 x 10“ 7 , and the resulting 
overall SER to approximately 2 x 10” 7 . The correspond- 
ing redundancy profile (60, 24, 34, 38, 40, 38, 34, 24) costs 
0.67 dB, which moves the concatenated code operating 
point to 0.87 dB. This is inferior to both the cases of one- 
in-four and the one-in- two strong codewords. 

All of the performance results discussed in this section 
for two-stage decoding have been obtained by trial-and- 
error optimization using Figs. 2 and 3, and they may not 
always reflect the exact optimum performance. However, 
because the concatenated code operating point stays rela- 
tively constant (within a couple hundredths of a dB) over 
a relatively large range of convolutional code operating 
points (at least a tenth of a dB) , the trial-and-error method 
gives results accurate to a few hundredths of a dB wherever 
data are available. 

C. Four-Stage Decoding 

Four-stage decoding consists of four passes each 
through the Viterbi decoder and the Reed-Solomon de- 
coder. The main difference relative to two-stage decoding 
is that now it is profitable to design four separate levels 
of Reed-Solomon codeword redundancies rather than two. 
The analysis proceeds as in the two-stage case, 

1. Stage VIT-1. The output of Stage VIT-1 is char- 
acterized by the performance curves given in Fig. 1. 

2. Stages RS-1, RS-2, RS-3, and and RS-4. The 
outputs of Stages RS-1, RS-2, RS-3, and RS-4 are charac- 
terized by the performance curves given in Fig. 2. 


3. Stages VIT-2, VIT-3, and VIT-4. The outputs 
of Stages VIT-2, VIT-3, and VIT-4 are characterized by 
the performance curves given in Fig. 3. 

4. Performance of Four-Stage Decoding. With 
four-stage decoding, it pays to design four different cat- 
egories of Reed-Solomon codewords, not just strong and 
weak. As an example, one codeword out of every eight 
could be assigned the highest redundancy, and these code- 
words could be expected to decode during Stage RS-1; one 
of eight could be assigned the next highest redundancy 
(spaced halfway between the highest redundancy words) 
and these could be expected to decode during Stage RS-2 
after Stage VIT-2 Viterbi redecoding, using a pattern of 
one known symbol every eight. Then the third highest 
redundancy words (one of every four spaced halfway be- 
tween the two types of higher redundancy words) will be 
decoded in Stage RS-3 after Stage VIT-3 Viterbi redecod- 
ing with the help of one known symbol every four; finally, 
the lowest redundancy words (one of every two) will be 
decoded in Stage RS-4 after Stage VIT-4 Viterbi redecod- 
ing using one known symbol every two. A scheme of this 
sort was proposed in [2] for the Galileo IIGA convolutional 
code and other codes. 

The overall SER for this scheme can be calculated as 
in the two-stage case by assuming that the SER for the 
codewords that are supposed to decode during a particu- 
lar stage also contributes to the SER for all weaker code- 
words. The resulting formula for the overall SER is of the 
form SER = SER a { 1) + 7/8 SER b { 2) + 3/4 SER C ( 3) + 
1/2 SERd( 4), where the indices a,6,c,d, refer to the 
strongest, next strongest, third strongest, and weakest 
codewords and (n) refers to decoding during Stage 
RS-n, n — 1,2, 3, 4. 

To achieve an overall SER of 2 x 10“ 7 , it is necessary 
that all four terms in the expression for SER be driven 
to approximately 10“ 7 or lower. At a convolutional code 
operating point of 0.10 dB, the equivalent operating points 
for Stages VIT-2 and VIT-3 are 0.33 dB and 0.72 dB, 
respectively, from the curves for phase ±4 of 8 and phase 
±2 of 4 in Fig. 3. An appropriate set of E values to drive 
SER a ( 1), SER h ( 2), and SER C { 3) each below 1 x 10“ 7 
is E = 36, 24, and 12. As in the earlier discussion, the 
SER-e quivalent output of Stage VIT-4 is beyond the scale 
of Fig. 2, but separate calcuations show that E — 4 is 
sufficient to bring SERd( 4) under 1 x 10“ 7 and, together 
with the three aforementioned values, to achieve an overall 
SER of about 2 x 10“ 7 . The corresponding redundancy 
profile (72, 8, 24, 8, 48, 8, 24, 8) costs 0.45 dB, and so the 
concatenated code operating point is E b /N 0 = 0.10 dB 4- 
0.45 dB = 0.55 dB. 
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If the convolutional code operating point is lowered to 

0.00 dB, the SER for phase ±1 of 2 is no longer small 
enough to permit E — 4 for the lowest redundancy code- 
words, but E — 5 is sufficient. From Fig. 3 the in- 
termediate SER - equivalent convolutional code operating 
points are 0.27 dB and 0.66 dB. An appropriate set of 
E values to achieve an overall SER of about 2 x 10“ 7 is 
E = 43, 27, 13, and 5. The corresponding redundancy 
profile (86, 10, 26, 10, 54, 10, 26, 10) costs 0.52 dB, 
and the concatenated code operating point is E^/Nq = 

0.00 dB + 0.52 dB = 0.52 dB. This represents a slight im- 
provement over the 0.10 dB convolutional code operating 
point, and there may be a slightly better operating point 
in the unsimulated region below 0.00 dB. 

IV. Analysis Based on Depth-8 Interleaving 
and No Erasure Declarations 

The analysis now continues for the actual Galileo condi- 
tions with depth-8 interleaving. As with the infinite inter- 
leaving analysis presented earlier, this section begins with 
a baseline analysis of one-stage decoding (no Viterbi rede- 
coding) and proceeds to evaluate the gains possible from 
using a second stage. Results are not presented for more 
than two decoding stages. 

A. One-Stage Decoding 

For one-stage decoding with depth-8 interleaving, the 
performance of the Viterbi decoder is the same as before. 
The performance of the Reed-Solomon decoder must be 
evaluated by directly measuring the probability of decoder 
failure on simulated interleaved output from the Viterbi 
decoder. 

1. Stage VIT-1. The output of Stage V1T-1 is char- 
acterized by the performance curves given in Fig. 1. 

2. Stage RS-1. Figure 4 shows the performance of 
the Reed-Solomon decoder under an assumption of depth- 
8 interleaving and no erasure declarations, with varying 
amounts of error correction E — 4, 10, 16, 32, and 34. The 
depth-8 interleaving curves are plotted against a backdrop 
“grid” of infinite interleaving curves from Fig. 2. 

In Fig. 4, the performance curves for depth-8 inter- 
leaving are shifted slightly to the right with respect to 
the corresponding infinite interleaving curves, showing a 
slight degradation due to nonideal interleaving of less than 
0.05 dB at SER values around 1 x 10~ 5 or so. Equiv- 
alently, this degradation may be characterized as an ap- 
proximate reduction in the error correction capacity of the 
Reed-Solomon code. For example, the depth-8 curve for 


E — 16 starts out coincident with the E = 16 infinite 
interleaving curve for high SER, but gradually drifts to 
cross the £7=15 infinite interleaving curve at an SER 
of - 1 x 10“ 3 , and appears headed to cross the E — 14 
infinite interleaving curve at an SER of ~ 1 x 10 -7 . Thus, 
depth-8 interleaving effectively decreases the error correc- 
tion capacity of the (255,223) Reed-Solomon code by ap- 
proximately two errors at an SER of ~ 1 x 10” 7 . 

For the curves in Fig. 4, as well as succeeding figures 
below, the results are only statistically meaningful down 
to about 1 x 10~ 5 or 1 x 10 -6 SER for the 2-Gbit datasets 
from Stage VIT-1. However, due to the almost parallel 
behavior of the depth-8 interleaving curves relative to the 
family of infinite interleaving curves, the depth-8 curves 
can be extended to 1 x 10 -7 SER with high confidence by 
extrapolating along the appropriate infinite interleaving 
curve(s). 

3. Performance of One- Stage Decoding. For the 
baseline coding system using the (255,223) Reed-Solomon 
code, E is 16 and a one-stage decoding system with depth- 
8 interleaving can achieve an SER of 2 x 10“ 7 (i.e., a 
BER of ~ 1 x 10~ 7 ) at convolutional code E h /N 0 = 
0.59 dB. This calculation assumes an extrapolation along 
the E = 14 infinite interleaving curve as discussed above. 
The corresponding concatenated code signal-to-noise ratio 
is E\>/N 0 = 1.17 dB, accounting for the 0.58 dB redun- 
dancy of the (255,223) Reed-Solomon code. 

B. Two-Stage Decoding 

For two-stage decoding with depth-8 interleaving, the 
performance of the Viterbi decoder in both stages is the 
same as that for the infinite interleaving analysis. The 
Reed-Solomon decoding performance is the same as that 
for one-stage decoding with depth-8 interleaving, to a de- 
gree of approximation noted below. 

1. Stage VIT-1. The output of Stage VIT-1 is char- 
acterized by the performance curves given in Fig. 1. 

2. Stage RS-1. The output of Stage RS-1 is charac- 
terized by the performance curves given in Fig. 4. 

3. Stage VIT-2. The output of Stage VIT-2 is charac- 
terized by the same performance curves as given in Fig. 3. 

4. Stage RS-2. The output of Stage RS-2 can be con- 
servatively approximated by the performance curves given 
in Fig. 4, after using the SER-e quivalent Stage VIT-1 op- 
erating points to characterize the output of Stage VIT-2. 
For the case of noninfinite interleaving, this equivalence 
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is not exact, because the output of Stage RS-2 depends 
not only on the input SER but also on the statistics of 
the error bursts that contribute to the average SER. The 
typical error bursts output from Stage VIT-2 were demon- 
strated to be much shorter than those output from Stage 
VIT-1 at the same value of SER. Shorter bursts are less 
likely to contribute two or more symbol errors to a given 
interleaved codeword, and hence the performance of Stage 
RS-2 should be slightly better (i.e., closer to infinite inter- 
leaving performance) than that of Stage RS-1 at the same 
average SER. 

5. Performance of Two- Stage Decoding. If one of 
every four codewords is strong, extrapolation of the curves 
in Fig. 4 shows that a correction capability approximately 
halfway between E = 32 and E = 34 for the strong words 
is sufficient to drive SER s (l) to about 1 x 10“ 7 at a con- 
volutional code operating point of 0.20 dB. According to 
Figs. 3 and 4, E — 10 is sufficient to push SER W ( 2) to 1 x 
10“ 7 for the weak codewords at phases ±1 and E = 11 
can achieve SER W { 2) = 2 x 10“ 7 for the weak codewords 
at phase ±2. The overall SER for this scheme is approxi- 
mately 1 x l(r 7 +(l/2)l x 10“ 7 + (l/4)2 x 10" 7 = 2x 10" 7 . 
The redundancy profile is (66,20,22,20) and costs 0.58 dB 
of overhead. The resulting concatenated code operating 
point is 0.20 dB + 0.58 dB = 0.78 dB. 


V. Analysis Based on Depth-8 Interleaving 
and a Simple Erasure Declaration Rule 

This section analyzes the gains possible from applying 
a combination of Viterbi redecoding and Reed-Solomon 
redecoding using erasure declarations to codewords inter- 
leaved to depth 8. The erasure declaration rule is a simple 
strategy modeled after one of the rules given by Belongie 
and Arnold 1 ; it is effective but not optimized. A symbol is 
always erased (a “double-sided” erasure) if it falls between 
two known erroneous symbols no more than eight symbols 
apart from each other. If this does not suffice to make the 
codeword decodable, a symbol is also erased (a “single- 
sided” erasure) if it is immediately adjacent to one known 
erroneous symbol. Single-sided erasures are extended one 
symbol farther (into the next adjacent codeword) every 
time another codeword successfully decodes. The erasure 
declaration process is repeated until no additional code- 
words are decodable within the block of eight interleaved 
words. 


1 M. Belongie and S. Arnold, “Error Forecasting Performance (Pre- 
liminary),” JPL Interoffice Memorandum 331-92.2-142 (internal 
document), Jet Propulsion Laboratory, Pasadena, California, Jan- 

uary 4, 1993. 


As with the analyses presented earlier, this section be- 
gins with a baseline analysis of one-stage decoding (no 
Viterbi redecoding) and proceeds to evaluate the gains 
possible from using a second stage. Results are not pre- 
sented for more than two decoding stages. 

A. One-Stage Decoding 

One-stage decoding consists of one pass through the 
Viterbi decoder followed by one or more passes through 
the Reed-Solomon decoder to allow the various sequences 
of erasure declarations described above to be tried. The 
performance of the Viterbi decoder is the same as before. 
The performance of the Reed-Solomon decoder must be 
evaluated by directly measuring the ultimate probability of 
decoder failure after applying the rule for erasure declara- 
tions, using simulated interleaved output from the Viterbi 
decoder. 

1. Stage VIT-1 The output of Stage VIT-1 is char- 
acterized by the performance curves given in Fig. 1. 

2. Stage RS-1. Figure 5 shows the performance of 
the Reed-Solomon decoder under an assumption of depth- 
8 interleaving and the simple erasure declaration algorithm 
described above. For the solid curves in Fig. 5, it is as- 
sumed that all eight codewords in a frame have the same 
error correction capacity E, and performance is plotted 
for E = 10, 16. The E = 16 performance curve in Fig. 5 
shows a performance gain of about 0.10 dB with respect 
to the corresponding curve in Fig. 4 and a gain of about 
0.05 dB with respect to the curve in Fig. 2. The erasure 
declaration algorithm is powerful enough not only to over- 
come the degradation due to finite interleaving but also to 
outperform decoding with ideal interleaving but no erasure 
declarations. At low SER , the E = 16 depth-8 interleav- 
ing curve with erasure declarations is almost coincident 
with the E = 18 infinite interleaving curve with no erasure 
declarations. Thus, depth-8 interleaving combined with 
Reed-Solomon redecoding using a simple erasure declara- 
tion rule effectively increases the error correction capacity 
of the (255,223) Reed-Solomon code at an SER of ~ 10“ 7 
by approximately two errors relative to infinite interleav- 
ing with no erasure declarations, or by approximately four 
errors relative to depth-8 interleaving with no erasure dec- 
larations. 

The dot-dash curve in Fig. 5 shows the performance of 
the Reed-Solomon decoder using the same erasure decla- 
ration algorithm but allowing the codewords to have the 
variable redundancy profile (44,28,28,28) discussed by Be- 
longie and Arnold. 2 This profile has the same average 

2 Ibid. 
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redundancy as the constant redundancy code with E — 16 
and therefore the same redundancy overhead. The SER 
curve for this variable redundancy scheme lies just under 
the infinite interleaving curve for E = 20. The effective 
average error correction capability is improved over the 
constant redundancy scheme by approximately two ad- 
ditional errors, or equivalently by about 0.05 dB at an 
SER of - 10- 2 * * 5 * . 

3. Performance of One- Stage Decoding. Extrap- 
olation of the E = 16 curve in Fig. 5 yields a required 
operating point of 0.46 dB to achieve an SER = 2x10 . 

This corresponds to a concatenated code operating point 
of 0.46 dB + 0.58 dB = 1.04 dB. The E = 10 curve re- 
quires 0.67 dB plus 0.35 dB of overhead or 1.02 dB for the 
concatenated code. The variable redundancy (44,28,28,28) 
code requires 0.40 dB + 0.58 dB = 0.98 dB. 

B. Two-Stage Decoding 

For two-stage decoding with depth-8 interleaving and 
Reed-Solomon erasure declarations, the performance of 
the Viterbi decoder in both stages is the same as be- 
fore. The evaluation of Reed-Solomon decoding perfor- 
mance must be modified somewhat to account for the fact 
that erasure declarations in strong codewords in Stage 
RS-1 must typically be forecast only from other strong 
codewords (because weak codewords are not highly likely 
to decode during Stage RS-1), and erasures in weak code- 
words in Stage RS-2 must be forecast from other weak 
codewords (because there are no symbol errors remaining 
in strong codewords in Stage RS-2 from which to forecast 
neighboring symbol errors). 

1. Stage VIT-1. The output of Stage VIT-1 is char- 
acterized by the performance curves given in Fig. 1. 

2. Stage RS-1. The output of Stage RS-1 is charac- 
terized by the performance curves given in Fig. 6. These 
curves show the SER that is obtained when strong and 
weak codewords are intermixed within a frame, but one at- 
tempts to decode only the strong words (and declares only 
those erasures based on decodable strong codewords). The 
plotted SER is the average SER obtained by accumulat- 
ing symbol errors from all the codewords in frames that 
contain at least one strong codeword error. The dashed 

curves in Fig. 6 show SER for a redundancy profile of 
one strong word every two codewords, and the thick solid 

curves show SER for one strong word every four code- 

words. The error correction values noted for these curves 

( E s = 26,28 for the first case and E s = 30,32,34 for the 

second case) pertain to the strong codewords only. 


3. Stage VIT-2. The output of Stage VIT-2 is char- 
acterized by the same performance curves given in Fig. 3. 

4. Stage RS-2. At present no curves directly char- 
acterize the output of Stage RS-2. The SER from Stage 
RS-2 can certainly be upper bounded by the SER ob- 
tainable without erasure declarations, as shown in Fig. 4. 
The incremental effect of Reed-Solomon redecoding with 
erasure declarations can be estimated to some degree by 
comparing curves in Fig. 4 (without erasure declarations) 
with corresponding curves in Fig. 5 (with erasure declara- 
tions), but this is not an exact procedure because the effect 
of erasure declarations is different after Viterbi redecoding 
than after Viterbi decoding the first time. In particular, 
the shorter bursts from Stage VIT-2 should lead to less ef- 
ficacious error forecasting than the longer bursts from the 
“equivalent” Stage VIT-1 operating point. 

5. Performance of Two-Stage Decoding. Assume 
that one of every four codewords is strong, and the convo- 
lutional code operating point is 0.20 dB. Extrapolation of 
the curves in Fig. 6 shows that a correction capability of 
E = 32 for the strong words is sufficient to drive SER s ( 1) 
to about 1 x 10 -7 . Similarly, according to Figs. 3 and 4, 
E = 10 is sufficient to push SER W ( 2) to the same level for 
the weak codewords at phases ±1 even without any assis- 
tance from erasure declarations. For the overall SER to 
reach 2 x 10~ 7 , the remaining weak codewords at phase ±2 
must achieve SER W {2) of 2 x 10' 7 or lower. From Figs. 3 
and 4 this can be approximately accomplished without any 
erasure declarations by a code with E = 11. However, it 
is estimated that this required value of E can also be re- 
duced to E = 10 by using erasure declarations from the 
codewords at phases ±1 to help decode the words in be- 
tween at phase i2. The redundancy profile for this case 
is (64,20,20,20) and it costs 0.56 dB of overhead. The 
resulting concatenated code operating point is 0.20 dB 
+ 0.56 dB = 0.76 dB. 

If the convolutional code operating point is moved 
to 0.10 dB, the corresponding required values of E for 
the strong and weak codewords are approximately E = 
38 and E - 11, respectively, but there are insufficient data 
plotted in Figs. 4 and 6 to confirm this directly. The corre- 
sponding redundancy profile is (76,22,22,22) and the con- 
catenated code operating point is 0.75 dB, which is almost 
identical to the previous value. Again this indicates an in- 
sensitivity of the optimum performance to small variations 
on the order of 0.1 dB in the convolutional code operating 
point. All other things being equal, there is a slight pref- 
erence to operate at the higher convolutional code E b /N 0 
in order to raise the channel symbol signal-to-noise ratio. 
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VI. Performance Comparisons 

Tables 1 and 2 summarize the optimal or near-optimal 
(within a few hundredths of a dB) concatenated code de- 
coding thresholds for achieving an overall SER of 2 x 10" 7 
(a BER of - 1 x 10 -7 ) for the various decoding alterna- 
tives considered in the preceding sections. 

Several conclusions can be drawn from these two ta- 
bles. Based on the infinite interleaving analysis, a second 
stage of Viterbi decoding without any Reed-Solomon re- 
decoding is worth about 0.38 dB at SER = 2 x 10“ 7 . The 
same comparison for depth-8 interleaving yields essentially 
the same answer, 0.39 dB. Adding Reed-Solomon rede- 
coding using erasure declarations increases this value only 
marginally to 0.41 dB for the case of depth-8 interleaving. 
This tiny incremental value of Reed-Solomon redecoding 
when used in conjunction with Viterbi redecoding is prob- 
ably somewhat underestimated due to the conservative ap- 
proximations made in the foregoing analysis, but the true 
incremental value is unlikely to be more than several hun- 
dredths of a dB. On the other hand, Reed-Solomon rede- 
coding without Viterbi redecoding is worth about 0.19 dB 
for depth-8 interleaving. 

Allowing four stages of Viterbi decoding (and four levels 
of Reed-Solomon redundancy) is worth about 0.58 dB for 
infinite interleaving and no erasure declarations. Data are 
not available to confirm the corresponding result for depth- 
8 interleaving, but a similar effect is expected. 


Tables 1 and 2 also give a measure of the amount by 
which depth-8 interleaving is nonideal. Compared to infi- 
nite interleaving, depth-8 interleaving costs 0.06 or 0.07 dB 
for both one- and two-stage decoding without erasure dec- 
larations. 


VII. Choosing a Coding System for Best Per- 
formance 

If the objective for the Galileo LG A code is to maximize 
the allowable data rate, this is accomplished by lowering 
the concatenated code’s decoding threshold as far as pos- 
sible. Under a constraint of no more than two decoding 
stages and two levels of codeword redundancy, the best 
performance is obtained by using two-stage decoding with 
one strong codeword every four. The redundancy profile 
(64,20,20,20) achieves a near-optimum concatenated code 
decoding threshold of 0.76 dB while requiring the convo- 
lutional code to operate at 0.20 dB. 

Based on the infinite interleaving analysis reported in 
Table 1, performance may be further improved by approx- 
imately 0.2 dB if four stages of Viterbi decoding and four 
levels of Reed-Solomon redundancy are permitted. Con- 
firmation of this effect and specification of the optimum 
four-level redundancy profile for depth-8 interleaving will 
be the subject of a future article. 
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SER AFTER REED-SOLOMON DECODING 






SER AFTER VITERBI REDECODING 



CONVOLUTIONAL CODE E b /N Q , dB 


Fig. 3. Symbol error rate after Viterbi redecoding with various 
spacings of known symbols and relative phases of unknown sym- 
bols. 


EQUIVALENT E^/Nq, FROM INITIAL VITERBI DECODING 


SER AFTER REED-SOLOMON DECODING 


CONVOLUTIONAL CODE E b /N 0 , dB 
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Fig. 4. Symbol error rate after Reed-Solomon decoding with depth-8 interleaving. 


109 




SER AFTER REED-SOLOMON DECODING AFTER REED-SOLOMON DECODING 


0.8 


0.9 1.0 


CONVOLUTIONAL CODE E b / N 0 , dB 



Fig. 5. Symbol error rate after Reed-Solomon decoding using erasure 
declarations with depth-8 interleaving. 


CONVOLUTIONAL CODE E b /N 0 , dB 



Fig. 6. Symbol error rate after Reed-Solomon decoding using erasure 
declarations among strong codewords only with depth-8 interleaving. 




Table 1. Performance comparisons for infinite interleaving (SE/7 =2X10 7 ). 


Decoding stages 

1 

2 

4 

Erasure declarations 

No 

No 

No 

Redundancy profile 

(32) 

(60, 18, 20, 18) 

(86, 10, 26, 10, 54, 10, 26, 10) 

Convolutional code E^/Nq (dB) 

0.52 

0.20 

0,00 

Concatenated code Et>/ No (dB) 

1.10 

0.72 

0.52 


Table 2. Performance comparisons for depth-8 interleaving ( SER — 2 X 10 7 ). 


Decoding stages 

1 

1 

2 

2 

Erasure declarations 

No 

Yes 

No 

Yes 

Redundancy profile 

(32) 

(44, 28, 28, 28) 

(66, 20, 22, 20) 

(64, 20, 20, 20) 

Convolutional code (dB) 

0.59 

0.40 

0.20 

0.20 

Concatenated code E^/Nq (dB) 

1.17 

0.98 

0.78 

0.76 


